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I review how hydrodynamical flow is related to the observed flow in ul- 
trarelativistic heavy ion collisions and how initial conditions, equation of 
state and freeze-out temperature affect flow in hydrodynamical models. 
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1. Introduction 

One of the goals of the experimental heavy ion program at ultrarelativis- 
tic energies is to create a dense system of strongly interacting particles. 
It is hoped that the particles formed in the primary collisions would 
rescatter often enough to reach local thermal equilibrium and behave as 
a particle fluid, not as a cloud of free particles. If such a state is reached, 
the finally observed particles should depict signs of collective behaviour 
such as flow. 

Our intuitive understanding of flow, i.e. collective motion, is closely 
tied to a classical macroscopic description of flow using the language 
and tools of hydrodynamics. This means that it is often easiest to use 
hydrodynamical concepts like temperature, pressure and flow velocity to 
describe collective motion even if the applicability of such concepts is far 
from certain. Hydrodynamical models are thus particularly suitable to 
describe flow phenomena, but we have to be careful not to confuse what is 
actually observed with our way of describing observations. An example 
of the limits of hydrodynamical language is that there is no generally 
accepted definition of flow in the context of heavy-ion collisions, but the 
word is used in its intuitive meaning. 

A rigorous definition of flow is beyond the scope of this review. In- 
stead I use a hydrodynamically practical definition of flow: collective 
flow is correlation of position and momentum during the dense, inter- 
acting stage of the collision regardless of the origin of these correlations. 
This means that I also call flow the correlation between the longitudi- 
nal momentum and the position of particles which has its origin in the 
initial particle producing processes. In a hydrodynamical model these 
correlations are manifested as an initial non-zero longitudinal velocity 
field. 

Unlike transverse and longitudinal flow, directed and elliptic flow do 
not directly refer to collective motion but to certain emission pattern 
where particle emission is not azimuthally isotropic (for definitions see 
sections l3.2.1l and l3.2.21 . In principle elliptic anisotropy could be entirely 
due to the shape of the surface of the source and be finite even if flow 
the velocity is zero^. Therefore to call elliptic anisotropy elliptic flow is 
unfortunate but firmly established in the literature. It has to be remem- 
bered, however, that even if elliptic anisotropy is not necessarily a sign 
of collective motion, it is a collective effect. 

Hydrodynamics connects the conservation laws to the equation of 
state, viscosity and heat conductivity of the fluid. Thus the properties 
of matter and flow are intimately connected and we hope to learn about 
the equation of state of nuclear matter by studying the flow in heavy ion 
collisions. In practice, however, this is a challenging task because of the 
nonlinear nature of the equations of hydrodynamics and the many un- 
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knowns in the hydrodynamical description of heavy-ion collision. In this 
review I describe briefly the basic concepts of a hydrodynamical model 
and the kind of collective flow generated in hydrodynamical simulation. 
My emphasis is on details which have particular significance in the de- 
scription of elliptic flow and how initial shape of the system, equation of 
state and freeze-out temperature affect elliptic flow in Au + Au collisions 
at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
Laboratory (BNL) at y/s = 130 GeV/A energy. 

The use of hydrodynamics has a long tradition in heavy ion physics. 
Cons equent ly there are many previous reviews and introductory ar- 
ticles ' ' ' ' where a reader can find a more detailed discussion of hydro- 
dynamical models. A complementary view of flow and its development in 
hydrodynamics is also provided in by Kolb and Heinz ^ in this volume. 

Flow in heavy-ion collisions is discussed from the experimental point of 

151 q 
view in reviews by Reisdorf and Ritter^and Herrmann et al.^ although 

their main emphasis is in collision energies below the y/s = 17.2 GeV/A 

collision energy of the Super Proton Synchrotron (SPS) at CERN. These 

reviews also provide a discussion of the experimental detection of flow 

which is beyond the scope of this review. 

2. Hydrodynamical model 

In a hydrodynamical description a heavy ion collision is basically as- 
sumed to proceed as follows: In the initial collision a large amount of the 
kinetic energy of the colliding nuclei is used to create a large number of 
secondary particles in a small volume. These particles will subsequently 
collide with each other sufficiently often to reach a state of local thermal 
equilibrium. When the system has reached local equilibrium it is charac- 
terized by the fields of temperature, T{x), chemical potentials associated 
with conserved charges, fj,i(x) and flow velocity, u^(x). The evolution of 
these fields is then determined by the hydrodynamical equations of mo- 
tion until the system is so dilute that the assumption of local thermal 
equilibrium breaks down and the particles begin to behave as free par- 
ticles instead. 

Besides the fact that we want to be able to describe the system using 
a few thermodynamic variables, the hydrodynamic description has ad- 
ditional advantages. Hydrodynamical models are relatively simple and 
one essentially does not need any information other than the equilibrium 
equation of state of nuclear matter to solve the equations of motion. Once 
the equation of state and the initial state of evolution are defined, the ex- 
pansion dynamics is determined and there is no need to know the details 
of the interaction on the microscopic level. This is especially practical 
when one wants to study the transition from hadronic to partonic degrees 
of freedom; a transition for which the details on the microscopic level are 



largely unknown. The use of familiar concepts like temperature, pressure 
and flow velocity also leads to an intuitive and transparent picture of the 
evolution of the collision. The price to be paid for these advantages is a 
set of bold assumptions: local kinetic and chemical equilibrium and lack 
of dissipation. This set of assumptions may or may not be valid in such 
a small system as that formed in a heavy ion collision. 

There is no proper proof for thermalization in heavy ion collisions. 
Instead, one has to revert to a heuristic comparison of collision rates of 
secondary particles with the lifetime of the collision system. At a temper- 
ature T>200 MeV, for example, the density of partons is n>4 fm -3 in a 
two flavor plasma. When the cross section is approximated by the pertur- 
bative QCD (pQCD) gluon-qluon scattering cross section of &gg~>gg ~ 3 
mb, the average mean free path, A = 1/an « 0.8 fm, and the time 
between two collisions is an order of magnitude smaller than hydrody- 
namically estimated lifetime of the system, r = 10-20 fm. Thus the 
partons should scatter several times and the system has a chance to 
thermalize. Another way to argue for thermalization in heavy ion colli- 
sions is simply to refer to the elliptic flow data which can be repro duc ed 
by a hydrodynamical model, but not as well by transport models^! 



2.1. Basics 

The equations of motion of relativistic fluid dynamics are the expres- 
sions for local conservation of energy and momentum and any conserved 
charge: 

d„T^ = and d M tf = 0, 

where T^ v is energy momentum tensor and , i = 1, . . . , n are the four- 
currents of conserved charges. Without any additional constraints these 
4 + n (n is the number of conserved currents) equations contain 10 + 4n 
unknown variables. The simplest and most commonly used approach to 
close this system of equations is the ideal fluid approximation which 
reduces the number of unknown variables to 5 + n. 

In the ideal fluid approximation the energy momentum tensor 

and currents ji are supposed to have forms 

T» v = (e + p)u^u v - pg^ and jf = , 

where e, p and n,; are energy density, pressure and number density of 
charge i in the local rest frame of the fluid and is the flow four- 
velocity of the fluid. In other words all dissipative effects such as viscosity 
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and heat conductivity are assumed to be zero and the fluid is always 
in perfect local kinetic equilibrium. The additional equation needed to 
close the system of equations is provided by the equilibrium equation of 
state (EoS) of the matter, which connects the pressure to the densities: 
P = P(e,ni,. . . , n n ). 

In principle it is possible to include small deviations from local ther- 
mal equilibrium by including dissipative effects, but in practice relativis- 
tic viscous hydrodynamics is very difficult to implement and h as not yet 
been done^l For estimates of the effects of viscosity, see refs. UAlSI 

The numerical solution of the hydrodynamical equations of motion 
in all three spatial dimensions is a tedious problem. In most approaches 
some approximate symmetry is applied to effectively reduce the number 
of spatial dimensions to two or one. A trivial simplification is to assume 
cylindrical symmetry in the description of central collisions of spherical 
nuclei. 

Another popular approximation is the Bjorken model l=™ where the 
longitudinal flow is assumed to be of the scaling form v z — z/t at all 
times. This requirement leads to boost invariance of the system: its pres- 
sure and energy density do not depend on the longitudinal coordinate 
z, if compared at the same proper time r = Vt 2 — z 2 . The solution of 
the equations of motion becomes independent of boosts along the beam 
axis and it is sufficient to solve the equations of motion in the trans- 
verse plane at z = 0. The obvious drawback in this approximation is 
that the results are independent of rapidity and one is limited to discuss 
on ly tran sverse behaviour. With the exception of results by Hirano et 
al. ' — "I the results discussed in this review are obtained in boost in- 
variant calculations. At RHIC energies the transverse flow results at mid 
rapidity are very similar in both boost invariant and non-boost invariant 
calculations. 



2.2. Initialization 

Local thermal equilibrium is one of the assumptions of a hydrodynamical 
model, but the model itself does not specify the mechanism that leads 
to an equilibrated state. Since at RHIC energies the initial particle pro- 
duction is definitely not an adiabatic process, hydrodynamics cannot be 
used to describe the initial collision, but the hydrodynamical evolution 
must begin at a sufficient time after the initial collision when the sys- 
tem has had time to reach thermal equilibrium. The initial state of the 
system, i.e. the density distributions and flow velocities at the beginning 
of the hydrodynamic evolution, are not given by the model either but 
must be given as external input. 

When analysing flow it is not enough to know what the maximum 
density or temperature reached in collision was. The maximum flow ve- 



locity before freeze-out depends also on pressure gradient close to the 
edge of the system. This makes particle distributions at high pr sensi- 
tive, not only to the maximum initial pressure, but also to the initial 
density profile ^1 The elliptic flow, on the oth er hand, is closely related 
to the initial asymmetry of the system G-7IJ-8 I This makes choosing the 
initial distributions an essential part of modelling flow. 

The simplest method to d eter mine the initial state is the one pro- 
posed by Hwa and Kajantie Since ideal fluid expansion is isen- 
tropic and entropy is related one-to-one to the measured multiplicity 
at a fixed freeze-out temperature and chemical potential, the final mul- 
tiplicity gives also the initial entropy of the system. Only a choice of the 
initial size of the system is needed to fix the (average) initial entropy 
density and if the equation of state is known, all other thermodynamical 
properties follow. This approach does not tell anything about the initial 
density distributions and more constraints are needed to study flow. 

In boost-invariant expansion it is sufficient to specify the density 
profiles in the transver se p lane. A plausible approach is to localize the 
wounded nucleon model'-*' and assume that the density in the transverse 
plane is proportional to the number of participants per unit area in the 
transverse plane. For two nuclei colliding with impact parameter b, the 
density of participants can be calculated from a geometric formula (ref.^ 
and references therein): 

Bl 



n WN (s;b) = T A (s+ ib) 



T B (s-ib) 



1 



VP 



T B (s- 



B 



a pp T A (s + j;b) 
A 



where a pp is the inelastic proton-proton cross section at the collision 
energy, T A is the nuclear thickness function of nucleus A; 

/oo 
dz p A {s,z), 
-oo 

and Pa{ f ) is nuclear density given by a Woods-Saxon distribution. 

At SPS, the final particle multiplicity scales with the number of 
participants ^4 Thus it is natural to assume that the initial entropy 
density scales with the number of participants J — L - iJ and the initial 
entropy density distribution is given by 

s(s;r ;b) = A' s (r )n WN (s; b), 

where K s (ro) is a proportionality constant chosen to reproduce the ob- 
served final particle multiplicity in central collisions. This constant de- 
pends on the initial time of hydrodynamical evolution, tq, which must be 
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chosen separately. In the following this parametrization is called sWN, 
as in ref . 

Strictly speaking the hydrodynamical approach does not allow both 
the transverse energy dEx/dy and the multiplicity dN/dy to scale with 
the number of participants. The data, however, does not rule out the 
possibility that transverse energy is proportion al to the n umber of par- 
ticipants'^. Thus one can make an assumption 12323211 that the initial 
energy density, not the initial entropy density, scales with the number of 
participants as given in Eq. 

e(s;r ;b) = A' £ (r )n WN (s; b), (7) 

with a proportionality constant K e ^ K s . By analogy to the label sWN, 
this initialization is labeled as eWN. 

With increasing collision energy one expects that the hard collisions 
between incoming partons become more and more important and finally 
dominate particle production 123 In that limit each nucleon-nucleon col- 
lision contributes equally to the particle and energy production and the 
number of produced particles scales like the number of binary nucleon- 
nucleon collisions in the transverse plane. It is given in terms of nuclear 
thickness functions (Eq. JSJ) as 

n BC (s;b) = o- pp T A (s + th)T B (s- ±b). (8) 

If the system thermalizes quickly via inelastic collisions the density of 
produced particles defines the initial entropy density at the beginning of 
the hydrodynamical expansion. Thus the initial entropy density should 
be proportional to the number of binary collisions: 

s(s;r ;b) = K s (To)n BC (s; b), (9) 

which defines parametrization sBC for initialization. The proportionality 
constant is now labeled k to emphasize that its value is different from 
the values of K s and K e in the sWN and eWN parametrizations. 

One can also argue that each binary collision contributes equally not 
only to particle production but also to the energy carried by the produced 
particles. In that case the initial energy density should be proportional 
to the number of binary collisions l 14 l 10 t 

e(s; r ; b) = K e (r )n BC (s; b). (10) 

This parametrization is called eBC. 

Each of these parametrizations leads to a different centrality depen- 
dence of the multiplicity. Thus one can use data to differentiate between 
them without any reference to flow. As shown in fig. the multiplicity 
data shows slightly stronger than linear dependence on the number of 



s 




Fig. 1. Charged particle yield per participating nucleon pair at midrapidity as a function 
of the number of participants for different initialization models discussed in the texil^l. 
All curves were normalized to dN„ f, /dn = 550 for 5% of the most central collisions 
(b = 2.3 fm.) The data are from refsPHEl 



participants I " I Parametrization eBC is closest to the data, but 

the linear behaviour of sWN is not far from the data either. Even if 
the multiplicity data constrains the initial parametrization somewhat, 
there is still freedom in choosing a combination of these parametriza- 
tions. One should also remember that these parametrizations are not 
the only possibilities. 

Besides the initial distributions one has to choose the initial time of 
hydrodynamic evolution, tq. As there is no method to calculate whether 
the system thermalizes, there is no method to calculate when the system 
is sufficiently thermalized for the hydrodynamic evolution to begin. Thus 
the initial time is another free parameter to be chosen to fit the data 
or by other arguments like saturation scale in pQCD calculations. In 
simulations of Au + Au collisions at RHI C (y^s = 130 GeV), the initial 
time has varied from tq = 0.2 to 1 fm™fl 

If the assumption of boost invariance is relaxed, the choice of initial 
state becomes considerably more complicated. There are few constraints 
for the longitudinal flow velocity profile or the longitudinal density distri- 
butions. Thus the choice of a particular parametrization and the values 
of the parameters is largely based on trial and error - tuning the model 
until a reasonable fit to experimental rapidity distributions is achieved. 
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For a sample of initial profiles used succe ssfully to describe longitu dinal 
expansion at the SPS see refs. I ^ l 33 l 34 l a nd at RHIC refs. IME4j It is 
also instructive to remember that even for the same EoS there are sev- 
eral pos sible initial states which lead to an acceptable reproduction of 
the data^. 

An alternative approach to determine the initial state is to use some 
other model to calculate it. For examp le event generator I^SI or pertur- 
bative QCD (pQCD) calculations GlEHl have been used for this purpose. 
Even if these approaches increase the predictive power of hydrodynam- 
ics, thermalization is still an additional assumption. 

2.3. Equation of State 

With the notable exception of ref. where the equation of state (EoS) 
is based on chiral SU(3) a — u model a all the hydrodynamical calcu- 
lations at RHIC energies have used equations of state based on similar 
structure: a hadronic phase which is constructed as a gas of free hadrons 
and resonances, a plasma phase of ideal, massless partons with a bag 
constant and a first order phase transition between these two phases. 

Different choices of the number of resonances included in the ha- 
dronic phase, of the latent heat and of the phase transition temperature 
cause minor differences in the equations of state of different practition- 
ers, but the major difference is whether the hadron phase is assumed 
to be in chemical equilibrium or not. Thermal models used to fit final 
state par ticle abundancies give larger freeze-out temperatures T c ^ ~ 160 
MeV than kinetic freeze-out temperatures Tf ~ 120 MeV obtained 
from fits to particle px spectra. This discrepancy can be explained by dif- 
ferent cross sections for elastic and inelastic collisions. Thermal models 
assume chemical equilibrium which requires frequent inelastic collisions 
which change particle number. On the other hand, frequent inelastic 
collisions are sufficient to maintain kinetic equilibrium. Since the cross 
sections for particle number changing collisions are much smaller than 
for collisions where particle number does not change (elastic and quasi- 
elastic collisions), it is natural to assume that inelastic collisions cease 
first and chemical freeze-out occurs at a higher temperature than kinetic 
freeze-out. Thus the system may be in local kinetic, but not chemical, 
equilibrium at the later stages of its evolution. 

Chemical non-equilibrium can be incorporated in the hydrodynam- 
ical description by treating lowest lying hadron states as stable parti- 
cles | 41 1 ^3 | rp ne p ar ^j c j e number of each of these hadrons forms a con- 



a In this paper only HBT radii were discussed and it is thus beyond the scope of this 
review. 
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served current and a finite chemical potential for each hadron is built 
up. Baryon and antibaryon chemical potentials are also independent in 
this approach. Chemical non-equilibrium changes the space-time evolu- 
tion of the system only slightly because the relation between pressure 
and energy den sity is very similar both in chemical equilibrium and non- 
equilibrium 1231. The main difference is that in chemical non-equilibrium 
the temperature decreases faster as energy density decreases and thus 
the system reaches its freeze-out temperature faster. How this changes 
the observed anisotropy will be discussed in section [I] 



2.4. Freeze-out 

At some point in the evolution the particles will begin to behave as free 
particles instead of a fluid and the hydrodynamical description must 
break down. When and where that happens is not given by hydrodynam- 
ics but must be included as external input. The conventional approach 
is to assume this to take place as a sudden transition from local thermal 
equilibrium to free streaming particles when the expansion rate of the 
system is larger than the collision rate between particles or the mean 
free path of the particles becomes larger than the system size. Finding 
where these conditions are fulfilled is a nontrivial problem. Since the 
scattering rate is strongly dependent on temperature the usual approxi- 
mation assumes that the freeze-out takes place on a hypersurface where 
temperature (or energy density) has a chosen freeze-out value. This tem- 
perature is of the order of the pion mass, but its exact value is largely 
a free parameter which can be chosen to fit the data. In Pb + Pb colli- 
sions at the SPS (y/s = 17 GeV/A) the calc ulat ed values of freeze-out 
temperatures vary between 100 and 140 MeV^I In Au + Au collisions 
at RHIC (^/s — 130 or 200 GeV/^4) the s pan of suggested freeze-out 
temperatures is even wider, from 100 ™ t0 160-165 MeV™. 

After choosing the surface where the freeze-out takes place, the ther- 
modynamic variables characterizing the state of the fluid must be con- 
verted to spectra of observable particles. A practical way of doing this 
is the Cooper- Frye algorithm "521 where the invariant momentum distri- 
bution of a hadron h is given by 

dN_ __ g h f 1 M 

* d*p (2^)3 J af expi(p^ - a)/T] ± 1 P ^ [ ' 

where the temperature T(x), chemical potential u(x) and flow velocity 
u^(x) are the corresponding values on the decoupling surface aj. Besides 
its relative simplicity, this approach has the advantage that if the same 
equation of state is used on both sides of decoupling surface, both energy 
and momentum are conserved. 
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However, the Cooper-Frye formula has a conceptual problem. For 
those areas where the freeze-out surface is spacelike, the product p^da^ 
may be either positive or negative, depending on the value and direction 
of p^. In other words, the number of particles freezing out on some parts 
of the freeze-out surface m ay b e negative. These negative contributions 
are small (a few per cent anc ] are usually ignored. More refined 
procedures without negative contributions have been suggested ^3 brrt 
their implementation is complicated. So far there has been no full-fledged 
calculation using these procedures. 

Another way to refine the hydrodynamical freeze-out procedure is 
to circumvent the entire problem and switch from a hydrodynamical to 
a microscopic transport model description well within the region where 
hydrodynamics is supposed to be applicable 1 — L - 21 . Besides giving a bet- 
ter description of freeze-out, such models include the separate chemical 
and kinetic freeze-outs. The main drawback of such models is - besides 
the increased complexity - that the correct region where the switch from 
the hydrodynamical to the transport description should take place is as 
uncertain as the kinetic freeze-out surface in ordinary hydro dynamical 
calculation. The educated guess employed in both refs.' — *— ^is that the 
switch happens immediately after hadronization. 



3. Transverse flow and its anisotropics 

In hydrodynamics, flow is generated by pressure gradients. In the initial 
state of a relativistic heavy-ion collision, there is only one type of col- 
lective motion, the coherent longitudinal motion of the two approaching 
nuclei. At center-of-mass energies above about 5GeV/A nucleon pair, 
this initial motion can no longer be completely stopped, even for large 
nuclei (^4 > 200) undergoing fully central collisions, and a certain frac- 
tion of the collective motion in the final state is simply a remnant of 
the initial beam motion. To separate it from hydrodynamically gener- 
ated longitudinal flow is notoriously difficult, and therefore most of the 
attention focuses on transverse collective flow in the directions perpen- 
dicular to the beam which was entirely absent before the collision and 
thus can be clearly associated with collective behaviour generated during 
the collision. 

One distinguishes between transverse flow in general and its anisot- 
ropics. In some works in literature, the azimuthally averaged transverse 
flow is called radial flow (e.g. ref. "J). However, in studies of heavy ion 
collisions at energies around 1 GeV / A, radial flow is understood to mean 
three dimensional spherically symmetric flow ^1^1 Therefore I prefer to 
call collective motion in the transverse plane transverse flow without any 
reference to its azimuthal structure. Transverse flow will be discussed in 
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Section 13.11 

In central (b — 0) collisions between spherical nuclei the transverse 
pressure gradients are independent of the azimuthal angle and transverse 
flow field is azimuthally isotropic. On the other hand in non-central 
(6 7^ 0) collisions between spherical nuclei, or in central collisions be- 
tween deformed nuclei (e.g. U+U), the nuclear overlap region is initially 
spatially deformed in the plane transverse to the beam, resulting in az- 
imuthal anisotropies of the pressure gradients and of the final collective 
flow pattern generated by them. Two specific forms of anisotropy of the 
transverse flow field are the "directed" and the "elliptic" "flows" dis- 
cussed in Sections 13.2.11 and 13.2.21 below. They are characterized and 
measured by the first two higher order expansion coefficients in an az- 
imuthal Fourier expansion of the final momentum spectra. These coef- 
ficients vary for different hadron species, but the hydrodynamic model 
relates their magnitudes in specific ways which depend on the initial and 
freeze-out conditions and the equation of state of the expanding matter. 

3.1. Transverse flow 

In one-fluid hydrodynamics, the concept of transverse flow is particularly 
simple. The pressure gradient between the dense center of the system and 
the ambient vacuum causes the system to expand and transverse velocity 
v r of the fluid is built up. However, if one defines the system size by the 
location of the freeze-out surface, the transverse flow velocity is not the 
expansion velocity of the system. Since the system dilutes, freeze-out sur- 
face may move inwards even if the particle fluid flows rapidly outwards. 
In such a case the whole concept of expansion velocity is ambiguous. 

The equation of state (EoS) is closely related to the buildup of flow: 
the stiffer the EoS, the larger the flow. Unfortunately many other factors 
affect the flow as well and therefore px spectra of particles constrain 
the EoS only weakly. Extrem e cases like the ideal pion gas EoS can 
be excluded as too stiff I^E BI but even the effect of a phase transition 
can be compensated by changes in the initial density and freeze-out 
temperature''"^. At large values of the transverse momentum, px > 2-3 
GeV/c, the particle distributions are also increasingly sensitive to the 
details of the velocity profile and thus to the initial pressure profile. 
Similar calculations may repr oduce the data or not depending on the 
choice of the initial profile'-*'. 

The experimental detection of transverse flow is much more difficult 
than its hydrodynamical description. Since it is not possible to detect 
where each particle was emitted, it is not possible to reconstruct the 
flow (or its absence) either. Instead one has to deduce the presence of 
flow indirectly by comparing p? distributions of various particle species. 
In the experimental literature the following procedure is often used to 
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Fig. 2. Slope parameters as a function of particle mass for (a) Pb+Pb central collisions 
at the SPS (<y/5 = ; 
GeV/A). From refi ; 



at the SPS (y/s = F7.2 GeV / A) and (b) Au+Au central collisions at RHIC (*fs = 130 
flM 



argue for collective flow: To characterize the slope of the px distribution 
the transverse mass spectra at midrapidity can be fitted to a simple 
Boltzmann distribution 

dN ( rn T \ 

oc exp . (12) 



Ay m T dm T \ T slopc 

Here r s i op c is the inverse slope parameter, often interpreted as the ap- 
parent temperature of the source. In p + p collisions at 450 GeV/cl^ 
the spectra of different particles (e.g. it, K, p) have a characteristic in- 
verse slope parameter of about 140 MeV. On the other hand, as shown 
in Fig. 1^1 in Pb + Pb collisions at the SPS ^ and in Au + Au collisions 
at RHIC the slope parameter of tt, K and p increases with parti- 
cle mass and collision energy. The linear increase in particle mass has 
been interpreted as a sign that T s i ope consists of two components: (a) 
thermal part, Tthcrm> associated with random motion, and (b) a part 
resembling collective motion with average transverse flow velocity ( v r )- 
Both contributions can be added and give rise to the slope parameter'" 

Tsiope = Tthcrm + m(v r } 2 . (13) 

However, one must warn against a too simpleminded use of this pro- 
cedure and interpretation of Eq. (I13H . First of all the experimental dis- 
tributions as well as the distributions given by Eqs. I|11I14|I are invariant 
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distributions whereas the thermal distribution used in the fits of Eq. 1121 
is not. The mj- factor required to make the distribution invariant and 
the integrations over the source leading to Bessel functions instead of 
exponentials cause the value of the slope parameter to depend on the px 
interval where the fit is carried out. 

To illustrate the emission from a boosted the rma l source we use the 
"blast-wave" model of Siemens and Rasmussen I™ where thermalized 
matter of temperature Tf, approximated by a boosted Boltzmann distri- 
bution, freezes out on a thin cylindrical shell. Assuming a boost-invariant 
longitudinal expansion, a transverse flow rapidity p on the shell and a 
freeze-out at constant proper time r, the Coope r-Frye freeze-out distri- 
bution (Eq. 1111 1 can be calculated analytically I 00 ! 00 ! Up to irrelevant 
constants one finds 

dN /p T sinhp\ /m T coshp\ . 

« m T I ^— — Z K x —!— , (14) 



dy mrp dmrp \ Tf J \ Tf 

where p = tanh -1 ty is the transverse flow rapidity. The slope parameter 



Tkope is given by^ 



T=^[^^-) d5) 



^siopc dm T \dym T dm T 

1 mj- sinh p cosh p 
m T PT T T 

where we have approximated the Bessel functions by exponentials. The 
slope parameter now has a much more complicated dependence on 
mass than the linear dependence of Eq. 1131 . As seen by differentiat- 
ing Eq. 1151 with respect to mass, 

dT s i p C _ 1 / m m sinhp\ 

dm ~ T2 \m^, m T p T T J ' [ ' 

there may even be a py-range where the slope parameter decreases as a 
function of mass. In a more realistic calculation the emission takes place 
on a surface where the flow velocity varies, and to find a px region where 
the slope parameter decreases with increasing mass requires a very large 
average flow velocity. At SPS and RHIC energies the flow velocity is 
not sufficiently large, but the slope parameters show a similar non-linear 
mass and pt dependence to that depicted in Eq. (1151 nevertheless. 

Besides the actual properties of boosted thermal distributions, there 
is another complication in extracting flow velocity from the px distribu- 
tions. A large number of detected particles do not originate from a ther- 
mal source but from resonance decays. Due to the available phase space, 
the decays contribute mainly to the low-py region leading to steeper 
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Fig. 3. Hydrodynamically calculated slope parameters as a function of particle mass in 
central collisions at RHIC (y/s = 130 GeV/A). In the left panel the slope parameters 
are from fits to spectra in two different p-r intervals. In the right panel the fits are done 
in interval 0.1 < Pt < 1 GeV/c for spectra after resonance decays (w. reso), before 
resonance decays (therm) and for thermal spectra immediately after hadronization at 
T c = 165 MeV. 



slopes than in the original thermal distribution. In ref. ^ the slope 
parameter of daughter particles originating in two body decays was ap- 
proximately related to the slope parameter of the resonance by 

Toft = — T R , (17) 
m R 

where p* , rag and Tr are the momentum of the daughter particle in 
the rest frame of the resonance, the mass of resonance and the slope 
parameter of the resonance, respectively. How much the decays change 
the observed slope parameter depends on the particular resonance and 
particle species and the temperature where the yields of resonances and 
particles are fixed. The calculation of the spectrum of resonance prod- 
ucts is generally a complex task and must be carried out numerically. A 
detailed discussion of decay kinematics is be yond the scop e of this article 
and an interested reader is referred to refs. I£Z1 58 | 59 | 60 | wnere different 
aspects of resonance decays have been discussed. 

How the fit intervals and resonance decays affect slope parameters 
can be seen in fig. [3] where hydrodynamically calculated slope parame- 
ters are shown. The calculation was tuned to reproduce the observed it 
and p spectra in the most central Au+Au collisions at \fs =130 GeV/A 
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energy at RHIC'-*'. The slope parameters in the left panel are from fits 
to the same spectra after resonance decays at freeze-out temperature 
at Tf — 130 MeV, but the px range where the fit was carried out was 
either 0.1 < pj- < 1 GeV or 0.5 < pt < 2 GeV. The pion and kaon slope 
parameters turned out to be quite independent of the fit interval, but 
all the heavier particles show a strong dependence on it. The slope pa- 
rameters in the right panel are all obtained using the former fit interval, 
but the spectra either contained the contribution from resonance decays 
(w. reso) or were the spectra of thermally emitted particles (therm). For 
comparison the slope parameters were also calculated for thermal spec- 
tra (Tj = 165, therm) immediately after hadronization at T c = 165 MeV. 
As explained, the resonance contribution makes the slopes steeper and 
decreases the slope parameter. This effect is small for heavy particles 
(m > irn) because there are very few resonances decaying into those 
particles. For the same reason, the mass dependence of the slope param- 
eter has a jump between the A and E masses. 

As seen in fig. |2]the experimental slope parameters of some particles 
do not follow the behaviour suggested by 7T, K and p. This has been 
used to argue that cj>, Q and J do not exp erien ce flow in the same 
way as pions and protons, but decouple earlier ^1 Since the scattering 
cross sections of 0-mesons and 57-baryons are smaller than pions and 
protons, this is possible, but as explained before, slope parameters are 
very sensitive to the vt interval of the fit and thus not very reliable. In 
their recent analysis the NA49 collaboration was able to fit the pt 
spectra of all particles using the blast wave model (Eq. ll4H and the same 
parameters for all particles including Q and Q. Thus the experimental 
situation is unclear as to whether all hadrons experience a similar flow 
or not. 

Nevertheless, the almost mass independent slope parameters shown 
in fig. |21 do not imply absence of flow. To illustrate this we show the 
slope parameters calculated from particle distributions immediately after 
hadronization in fig. [3] Even if there is flow at that time, the slope 
parameters show only a weak dependence on mass because the average 
flow velocity is small (0.2 vs. 0.44 at the end of hadronic phase) and the 
slope parameters are dominated by the large temperature. 



3.2. Flow anisotropics 

The anisotropy of transverse flow is manifested as azimuthally anisotro- 
pic final particle distribution. To quantify this anisotropy, the particle 
spect ra a re expanded in harmonics of the azimuthal angle cfi event by 
event ^ 

dN dN 



dyd(j)p 2-Kdy 



(1 + 2vi cos((/> — <f>ji) + 2v2 cos 2(0 — 4>r) 
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Fig. 4. Schematic representation of the collision geometry and different anisotropics 
of flow seen in the transverse plane. P and T denote the projectile and target nuclei, 
respectively. Top: directed flow in projectile rapidity region, positive (left) and negative 
(right). Bottom: elliptic flow, in plane (left) and out of plane (right). From ref" 



diV 



dN 



dydpT<i(j>p 2irdydpT 



(1 + 2vi(p T ) cos(eS - <f> R ) 
+2v 2 (pt) cos 2(<f> - (f> R ) + ■ 



(18) 



where (j>fj is the azimuthal angle of the reaction plane. Assuming that 
the experimental uncertainties in event plane reconstruction can be cor- 
rected for, each event can be rotated such that 4>r = 0. The expansion 
parameters t>i and v 2 correspond to directed and elliptic flow, respec- 
tively. Due to symmetry both vanish in central collisions which are cylin- 
drically symmetric. In symmetric collision systems the odd coefficients 
vi, V3, . . . , are zero at midrapidity. 



3.2.1. Directed flow 

At AGS energies the general picture of the origin of directed flow ^ is 
that the pressure formed in the collision region deflects the projectile and 
intermediate rapidity fragments, i.e. spectator nucleons away from the 
target ("bounce-off" and "sidesplash" effects^!) resulting in a preferred 
direction of nucleon emission. Simultaneously the produced pions scatter 
from spectator nucleons forming resonances. Through this process the 
initial direction of pions is lost and flow into the direction of spectator 
matter is reduced. Thus pions show flow to the opposite direction than 
nucleons. The usual sign convention is to choose protons to have positive 
directed flow in the projectile rapidity region and negative in the target 
rapidity region. Directed flow is zero at midrapidity due to symmetry 
and saturates before projectile and target rapidities are reached forming 
an overall s-shaped curve as function of rapidity^. 

The strength of directed flow depends on the pressure formed in the 
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collision region and the time that the secondary particles have to inter- 
act with the spectators, i.e. how fast the spectators pass by the collision 
zone. This time gets smaller with increasing energy and one expects the 
directed flow at SPS be smaller than at AGS, as is observed^. The for- 
mation of directed flow in the very early stages of the collision in very 
short times poses a problem for the hydrodynamical calculation. Espe- 
cially at SPS and RHIC energies it is possible that most of directed flow 
is established before the system has reach ed local thermal equilibrium 
and the pre-equilibrium features dominate ^! 

So far there are no detailed hydrodynamic calculations of directed 
flow at SPS or RHIC energies. However, it has been predicted that a 
phase transition in a nuclear equation of state would lead to a local 
minimum of directed flow as a function of collision energy "I The exact 
position and value of this minimum depend on the details of the model, 
but it should be located between the AGS and the maximum SPS energy. 

Another interesting prediction is that when the collision energy in- 
creases, the rapidity dependence of directed flow would eventually have 
such a shape that the si gn of direc ted flo w ch anges three times as a 
function of rapidity T n re f s . H3 and^ this behaviour is called 

"antiflow" and "third flow component", respectively, and its origin is 
explained by the phase transition and the initial shape of the system. 
The phase transition leads to an initial shape of a disk slightly tilted 
away from the beam axis. Thus the largest pressure gradient and the 
direction of fastest expansion points away from the spectator nucleons 
and the direction of the conventional directed flow. Close to midrapidity 
the emission from the collision region dominates and the particles show 
negative directed flow in the forward and positive in the backward rapid- 
ity region, whereas close to fragmentation regions directed flow is again 
due to interactions wit h sp ectator matter and has its usual sign. The 
initial ar gum ent of ref. I™ was for AGS energies, but has been refined 
for RHICI23 

Similar behaviour for a completely different reason is predicted in 
ref. "I where the changes of sign are predicted to happen above SPS 
energies independent of a phase transition. Their argument is that in 
non-central collisions nuclear transparency leads to an anisotropic dis- 
tribution of baryon number in the initial state, which is then manifested 
as anisotropic particle emission: when the particle fluid expands, baryons 
are carried to the same side as they were at the beginning of expansion. 
Even if the authors do not mention it, these two approaches can be 
expected to have obse rvable differences in pion directed flow: emission 
from a tilted disk ^01 i s similar for both pions and protons and both 
should show directed flow in the same direction close to midrapidity and 
in opposite directions in fragmentation regions. On the other hand, an 



Hydrodynamical description of collective flow 



inhomogeneous source'*^' leads to pion flow in the opposite direction to 
proton flow at all rapidities, or alternatively to a negligibly small directed 
flow of pions. From the hydrodynamical point of view, both of these ap- 
proaches are possible which highlights the uncertainties in choosing the 
initial state of hydrodynamic evo lution (section 12. 2|l . 

Unfortunately the SPS data c\ oes no t have a large enough 

rapidity coverage to test these ideas. So far there is no measurement of 
directed flow at RHIC. Whether it is large enough to be experimentally 
observable remains to be seen. 



3.2.2. Elliptic flow 

In non-central collisions the primary particle production is azimuthally 
isotropic, but since the interaction region is anisotropic in space, the 
secondary collisions can cause the final particle distribution to be anisot- 
ropic in momentum space. As depicted in fig.|S] particles heading out of 
plane have on the average a longer distance to go within the dense region 
than particles moving in plane. Thus particles moving out of plane have 
a larger probability to scatter several times and change their direction 
than particles heading in plane. Thus the final particle distribution has 
more particles moving in the in plane than the out of plane direction 
and the Fourier coefficient V2 in Eg. 1181 is positive. 

In a hydrodynamic picture the buildup of such elliptic anisotropy 
can be understood in terms of pressure gradients. The average pressure 
gradient between the center of the system and the surrounding vacuum is 
larger in plane than out of plane direction because the system is thinner 
in that direction. Consequently the collective flow velocity is larger in 
plane than out of plane. This leads to a larger average momentum in 
plane than out of plane and to more particles being emitted in the in 
plane than the out of plane direction. 

The origin of elliptic anisotropy is thus rescatterings and the shape of 
the system. The asymmetry of the shape of the system is largest immedi- 
ately after the collision and decreases with increasing time independently 
of the frequency of rescatterings. Only the speed with which the system 
will gain an azimuthally symmetric shape depends on how frequent the 
rescatterings are. Thus a large elliptic anisotropy and a large value of 
V2 is taken to be a signal of abundant rescatterings in the early stage of 
the collision and thus a signal of early pressure buildup and thermaliza- 
tion I2SIS3 Since hydrodynamics assumes zero mean free path and thus 
infinite rescattering, it is also assumed to give the practical upper limit 
of elliptic anisotropy. 

The relation between large vi and early pressure buildup is difficult 
to quantify because of the uncertainty of the initial time of the hydrody- 
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namical evolution. If pressure buildup and thermalization take a "long" 
time, say several fm/c, the particles formed in the primary collisions 
can move significantly in the transverse plane. The geometric arguments 
presented in section IT^l to constrain the initial shape are no longer valid 
and the initial state of hydrodynamical evolution is largely unknown. 
Consequently it is not possible to calculate reliably the elliptic flow pa- 
rameter V2 as function of thermalization time tq. The uncertainty of 
shape also explains the conflicting results in th e lit erature: if the shape 
does not change, vi is almost independent of tq^I On the other hand, 
if the shape changes as if the particles were freely streaming and the 
ratio of the anisotropy of the initial shape and the final V2 stays the 
same, a del ay of 3.5 fm/c in thermalization was estimated to reduce v<i 
by 50% EH. 

To gain insight how elliptic anisotropy is built up during the hydro- 
dynamical evolution one uses so called momentum anisotropy!^, labeled 
either 8 or e x in the literature, 

_ {T xx ~ T yy ) _ J dx dy (T xx - T yv ) 
€x ~ ~ (r xx + Tvv) ~ J dx dy (T xx + Tw) ' 

where T xx and T yy are the components of energy-momentum tensor 
T A " y , the angle brackets denote averaging over the transverse plane and 
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Fig. 6. Time-evolution of momentum anisotropy (solid line), elliptic anisotropy of ther- 
mal pions (dashed line) and a sum of thermal pions and pions originating from resonance 
decays (dotted line) in Au+Au collision at y/s = 130 GeV/A with impact parameter 
b = 7 fm. 



integration is done over transverse plane at constant proper time r. As 
an example the time evolution of momentum anisotropy (solid line) in 
a simulated Au+Au collision at \fs = 130 GeV/A collision energy is 
shown in fig. |S| 

The shoulder in th e inc rease of momentum anisotropy can be related 

to the phase transition^! In this particular case it occurs when no part 

of the system is in the plasma phase anymore, but the anisotropy begins 

|24l 

to increase before the system is entirely hadronized. As shown in ref. 1 — 
the actual effect of a phase transition on the buildup of momentum 
anisotropy is complicated and depends on the details of the flow field 
at the time of the phase transition and the relative sizes of the plasma, 
mixed and hadronic phases at each point of time. It is possible that the 
increase of anisotropy is halted as in fig.|Hl but it is equally possible that 
anisotropy keeps increasing or even decreases during the mixed phase. 
The decrease during the mixed phase can be explained by noting that 
the pressure gradients within the mixed phase are tiny and therefore 
the matter in the mixed phase keeps flowing with the flow velocity it 
had when entering the mixed phase. When the matter flows outwards 
with constant velocity, its flow is nearly self-similar, which decreases the 



momentum anisotropy. Whether this effect can be seen in the momentum 
anisotropy of the entire system, depends on how large a part of the 
system is in mixed phase. 

Another feature of the time evolution of the momentum anisotropy is 
that after saturating at t — to ~ 10 fm, it begins to slightly decrease. This 
is a result of the self-cancelling nature of the anisotropy; the anisotropic 
pressure gradients cause the system to expand more strongly in-plane 
than out of plane. In the course of time the system shape becomes cylin- 
drically isotropic, but due to momentum conservation the larger expan- 
sion in-plane continues. The system becomes elongated in-plane and the 
anisotropy of the pressure gradient begins to act against the momentum 
anisotropy. 

The relationship between the momentum anisotropy, e p , and the ob- 
served elliptic anisotropy, v%, is nontrivial. The former measures the 
anisotropy of the collective flow velocity, whereas the latter the anisot- 
ropy of particle yield. The particle distributions are a result of both col- 
lective and thermal motion and therefore the ratio of e p and vi depends 
on the mass of the particle in question and the freeze-out temperature. 
To illustrate this dependence we show in fig. |S| also «2 for pions calcu- 
lated assuming freeze-out at various temperatures 165 < Tf < 90 MeV 
and plotted as function of the freeze-out time of the system (dashed 
line). It can be seen that the momentum anisotropy increases roughly 
by 25% after the system is completely hadronized (at r — tq « 6.5 fm/c), 
whereas the anisotropy of pions increases by 70% - the difference caused 
by the increasing sensitivity of particle distributions to collective motion 
when temperature decreases. 

Another complication in comparing the momentum anisotropy of 
the flow to the observed anisotropy of particle distributions are the res- 
onance decays. The kinematics of decay not only favours the daughter 
particles having a smaller pj- than the decaying resonance, but also an 
azimuthal distribution of daughter particles which is peaked out-of-plane 
even if the distribution of the resonances peaks in-plane^! This leads to 
daughter particles showing negative elliptic anisotropy, which dilutes the 
anisotropy of particles with thermal origin, as can be seen in fig. H3 where 
the dotted curve shows the vi of pions originating both from decays and 
a thermal source. The effect of resonance decays depends on tempera- 
ture because the relative abundances of particles and resonances depend 
on temperature - the higher the temperature, the larger the fraction of 
pions that originates from resonance decays. 

We also show the ratio e p fvi as function of temperature in fig. Q 
It is seen that vi for thermal pions is roughly half of the momentum 
anisotropy, but when resonance decays are included, the ratio is much 
more temperature dependent. 

In general the heavier the particle, the more sensitive it is to collec- 
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Fig. 7. The ratio of momentum anisotropy, e p , to elliptic anisotropy, V2, of pions as 
function of freeze out temperature Tf in Au+Au collisions at ^/s = 130 GeV/A with 
impact parameter 6 = 7 fm. 

tive motion. For example, as explained in chaptcr l3.il if the freeze-out 
temperature and flow velocity are the same, the heavy particles have flat- 
ter px distributions and thus larger slope parameters T s \ ope than lighter 
particles. In the same way, the px averaged elliptic anisotropy parameter 
V2 is observed to increase with particle mass. However, at low values of 
PT, the pT-differential anisotropy parameter «2(pt)> shows the opposite 
behaviour: the heavier the particle, the smaller the value of V2(pt) at 
fixed px (see fig.[8Jl. 

The apparent contradiction between V2 and i>2 (pt) nas a simple ex- 
planation. V2 is not an additive quantity, but when a p^-averaged V2 is 
calculated from V2(px)y the latter is weighted by the particle distribu- 
tion: 



Thus the flatter px distribution of heavier particles weights the high px 
region, where V2(px) is larger, more and the p^-averaged value can be 
larger. Whether this larger weight for high px wins over the reduction 
of V2 at fixed px depends on the details of the expansion dynamics and 
the contribution from resonance decays. 

The mass dependence of the elliptic anisotropy at fixed px can be ex- 
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plained as an interplay between transverse co llec tive flow, random ther- 
mal motion and anisotropy of the flow field ^1 It is well known that 
transverse flow shifts the p^-distributions to larger values of px- For 
nonrelativistic px < m this effect increases with the particle mass m 
and the transverse flow velocity (v±). In the extreme case of a thin shell 
expanding at high velocity, the spectrum actually develops a rel ativ e 
minimum at px = and a peak at nonzero px ("blast wave peak" 
and with increasing mass the peak shifts to larger px- Relative to the 
case without transverse flow, the spectrum is thus depleted at small px, 
and the depletion, as well as the px range over which it occurs, increase 
with m and 

If the transverse velocity is larger in the x than the y direction, the 
same is true for this relative depletion effect. It counteracts the overall 
excess of particles moving to the x direction over particles moving to the 
y direction reducing V2- This reduction and the range where it occurs 
increases with particle mass and average transverse flow velocity In 
the extreme case of a thin shell, this depletion effect can be so large that 
there are less small px particles moving in the x than the y direction, 
and V2 becomes negative at low px ■ When a constant expansion velocity 
of a thin shell is replaced by a realistic transverse velo city distribution, 
the peak in the single particle spectrum disappears^! In the same way 
a realistic velocity profile weakens the reduction of vi at low px, but 
the mass dependence of V2 at low px remains. Whether particles show 
a positive or negative vi at low px depends on the details of the flow 
profile. 

For relativistic px > m, the particle mass does not play a role in 
the thermal distribution, and consequently the curves showing «2(pt) 
approach each other. In a simple model where the transverse velocity 
profile is replaced by its av erage value, as in theblast wave model of 
Siemens and Rasmussen 1^ and its derivatives l ^6 | 43 ] v ^ mcreases w ;th 
Px and approaches an asymptotic value of one. The details of the velocity 
profile can change this behaviour, but so far hydrodynamical calculations 
with realistic initial conditions have shown similar monotonic increases 
of V2 with px (see fig. [SJ . 



4. Learning from RHIC data 

One of the first measurements of Au + Au reactions at the Relativistic 
Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) 
at ^J~s — 130 GeV / A energy was a measurement of elliptic flow at midra- 
pidity Soon afterwards the elliptic flow in minimum bias collisi ons 
was shown to be very nicely reproduced by a hydrodynamical model ^1 
As shown in fig. |H] the differential anisotropy V2(px) m minimum bias 
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P t (GeV/c) Pt (GeV/c) 

Fig. 8. The calculated pt differential elliptic anisotropy coefficient V2(pr) at midra- 
pidity in minimum bias Au+Au coll isions at t/s = 130 GeV/A. The left figure shows 
v 2(pt) f° r various particle spec ie J^SI and the ri ght figure a comparison of the calculated 
V2{pt) for charged particleJlSl with STAR datsESl 



collisions follows the hydrodynamical curve closely up to momentum 
PT = 1-1.5 GeV As a function of centrality the hydrodynamical descrip- 
tion works as well or even better up to 16% of the most central collisions 
(6<6 fm), see fig. El At larger impact parameters, the px region where 
the hydrodynamical calculation fits the data becomes smaller and the 
deviation from the hydrodynamical curve grows faster. The identified 
particle anisotropy shows a similar mass ordering to that discussed at 
the end of previous section (figs. llOl and lllH . although the quantitative 
fit, especially of kaons, is not as good as the fit of pions and charged 
particles. The system thus behaves as a thermal sys tem and this result 
has been taken as an indicator of thermalization 

Since flow is closely connected to the equation of state of matter, 
it is possible to use details of the flow to constrain it. However, flow is 
equally sensitive to freeze-out temperature and the initial state of the 
evolution, so to learn about equation of state one must also know how 
flow depends on all other variables. 

The basic rule as to how px averaged anisotropy depends on the 
freeze-out temperature at RHIC energies is simple. The smaller the 
freeze-out temperature, the larger the i>2 — if everything else in the 
simulation stays unchanged. This rule is not valid indefinitely. As ex- 
plained in Section 13.2.21 and seen fig. |S] when the shape of the system 
has become azimuthally symmetric, the pressure gradients begin to work 
against the buildup of anisotropy and the momentum anisotropy begins 
to decrease. The observed i»2 may still keep increasing after the mo- 
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mentum anisotropy has begun to decrease because lower temperature 
makes the final particle distributions more sensitive to the anisotropy of 
collective motion, but eventually it will begin to decrease too. 

Unfortunately the temperature dependence of the px differential 
anisotropy, U2(pr)i is much more complicated. It is different for each 
particle species and it is also sensitive to the equation of state used. 
As shown in Fig. 1111 a decrease of freeze-out temperature from 140 to 
120 MeV, increases the pion V2(pt) slightly! but decreases the proton 
v 2(pt)- This again is a result of an interplay of collective flow and ran- 
dom thermal motion which is different for particles of different mass. 
The change in freeze-out temperature means two things: first, the life- 
time of the system changes and therefore the amount of flow changes. 
Second the random thermal motion changes. Since the observed particle 
distributions are the result of both, the changes in V2 result from both. 
Intuitively the effect of freeze-out temperature can be explained using 
the blast wave model in the same way as the mass dependence of v% (px) 
(Section l3.2.2H 'l. When freeze-out temperature decreases, transverse flow 
velocity increases and the blast wave peak in the proton px distribution 
becomes more prominent and moves to larger px- This makes the rela- 
tive depletion effect at low px larger and thus decreases the anisotropy. 
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Fig. 10. The differential elliptic anisotropy i>2(pt) °f kaons and Lambdas in minimum 
bias Au + Au collisions at ^/s = 130 GeV/A compared with the STAR dat i 43 ! 78 ! The 
calculation is similar to that of ref. 



The pion px distribution on the other hand does not show a similar peak 
or it lies at very small values of pt- Therefore there is no such depletion 
effect which would decrease «2 or it is limited to very low values of px- 
Fig. Illl also shows «2(pt) calculated using equations of state (EoS) 
with a phase transition (Q) and without a phase transition (H). In both 
cases the change in freeze-out temperature changes i>2(pt) as described. 
This is not surprising since below the phase transition temperature T c = 
165 MeV of EoS Q, both equations of state are identical. On the other 
hand, if the hadronic part of the equation of state is different, a change in 
the freeze-out temperature may affect «2(pt) differently. This is shown in 
fig. 1121 where V2(pt) f° r pions, kaons and protons is calculated assuming 
either local chemical equilibrium until kinetic freeze-out (CE), or that the 
relative abundancies of each particle species was frozen out immediately 
after phase transition at T c ^ = 170 MeV (PCE). In other words, the 
chemical freeze-out takes place at much larger temperature than kinetic 
freeze-out. 

In the case of chemical equilibrium (CE) the equation of state is 
very similar to EoS Q used to calculate the results shown in fig. Illl and 
the freeze-out temperature dependence is similar. The only difference is 
that when temperature decreases from 120 MeV to 100 MeV, the pion 
v 2(pt) does not increase anymore but begins to decrease. If chemical 
equilibrium is lost already at hadronization, the behaviour changes. The 
proton i)2 (pt) stm decreases with temperature, but the dependence is 
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Fig. 11. The differential elliptic anisotropy VzIpt) of pions and protons in mini mum 
bias Au + Au collisions at *Js = 130 GeV/iJ^l compared with the STAR dat;E3 The 
letters Q and H in the labels stand for an equation of state with a first order phase 
transition and a hadron gas equation of state without a phase transition, respectively. 
Numbers in parentheses stand for the freeze-out temperature in MeV. 

much weaker. On the other hand the pion i>2(pt) increases when Tj 
decreases in the entire temperature interval 100 < Tj < 140 MeV and 
the increase is much stronger than in the case of chemical equilibrium. 
The most dramatic change is in the kaon «2(pt) : if chemical equilibrium 
is maintained, i>2(pt) decreases when Tj decreases, but in the case of 
non-equilibrium, the behaviour is opposite. The kaon V2(pt) behaves 
like that of pions and increases as temperature decreases. Thus there 
are no simple rules to tell how the observed V2{pt) would correspond to 
a certain freeze-out temperature, but the constraints must be searched 
case by case for each equation of state. 

In general a stiffer equation of state causes larger flow. This holds also 
for elliptic anisotropy and, at least in the case shown in fig. 1111 also for px 
differential anisotropy: EoS H leads to a larger «2(pt) f° r both pions and 
protons than EoS Q. For the chemical equilibrium or out-of-equilibrium 
equations of state in fig. H2l this rule is less clear because both equations 
of state are almost equally stiff — the pressure as a function of energy 
density is almost identical in both cases. The difference between CE and 
PCE comes from the relation between temperature and energy density. 
The same energy density corresponds to a smaller temperature in PCE 
than in CE and therefore the system cools faster. Correspondingly the 
flow at fixed temperature is smaller when the system is out of chemical 
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Fig. 12. The differential elliptic anisotropy t>2(pr) of pions, kaons and protons in 
minimum bias \ Au + Au collisions at yfs = 13 G eV/A at different kinetic freeze-out 
temperature^^! compared with the STAR datJSl The left and right panels represent 
the results calculated assuming local chemical equilibrium until kinetic freeze-out (CE) 
or chemical freeze-out already at T c h = 170 MeV (PCE). 



equilibrium and the final anisotropy looks very different in these two 
cases. Change to chemical non-equilibrium means ending the evolution 
at an earlier stage, but the effect is not as straightforward as increasing 
the freeze-out temperature, because the random thermal motion is not 
changed. 

Constraining the equation of state by the transverse momentum 
spectra alone is notoriously difficult since larger flow generated by a 
stiffer equation of state can to a large extent be compensated b y a 
higher freeze-out temperature and a slightly different initial state 
As can be seen in fig. 1111 essentially the same holds for pion elliptic flow: 
a purely hadronic equation of state (H) and an equation of state with 
a phase transition (Q) create similar «2(pt) f° r pions if the freeze-out 
temperature is chosen to be Tj = 140 and 120 MeV for EoS H and Q, re- 
spectively. On the other hand, the effect on the proton i>2(pt) ' s exactly 
opposite: this choice leads to the largest difference of all the combina- 
tions studied. The different sensitivity of the proton and pion i>2(pt) 
to the equation of state and freeze-out temperature gives an additional 
handle on constraining both. 

The results in fig. llll cleaxlv favour an equation of state with a phase 
transition. As well one could claim, based on fig. 1121 that the proton 
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Fig. 13. The differential elliptic anisotropy V2(pt) of p ions and protons in minimum 
bias Au + Au collisions for different initialization modelJi^l compared with the STAR 
data 4 3. 



v 2(pt) favours an equation of state where chemical equilibrium is main- 
tained for temperatures lower than to T c ^ = 170 MeV. Unfortunately 
both conclusions are premature since there are other variables to be 
considered which affect the observed differential anisotropy. 

As explained in Section l2"l2l there are various ways to parametrize the 
initial state of hydrodynamical evolution. The simplest ones described 
in Section [2.21 all lead to slightly different shapes of the initial system. 
If the freeze-out temperature is kept unchanged, the ratio of the initial 
spatial anisotropy and the observed elliptic anisotropy is almost inde- 
pendent of impact parameter in a hydrodynamical model E' l ^^ l Despite 
that, the initial shape has its effect on the differential anisotropy ^(pt) 
in minimum bias collisions. The anisotropy calculated using the in itial 
state parametrizations described in Section [2.21 are shown in fig. [SI 18 . 
The pion anisotropy shown in the left panel is practically independent 
of initialization, but the proton V2(pt) is clearly sensitive to the initial 
shape. These results shown in fig. 1111 were obtained using eWN initial- 
ization and in principle it is possible that combining sBC initialization 
with EoS H would the reduce proton anisotropy sufficiently to reach the 
data. 

The most comprehensi ve s tudy of pt spectra and anisotropics at 
both SPS and RHIC so far reached the conclusion that an equation 
of state with a phase transition is necessary for a consistent reproduc- 
tion of the data. This study used a hybrid model of hydro and transport 
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Fig. 14. Pscudorapidity dependence of elliptic anisotropy in mini mum bias Au+Au 
collisions at s/s = 130 GeV/A compared to STAfEHI and PHOBO^HSl data assuming 
chemical freeze out at kinetic freeze out (CE) or at hadronization (PCE)I3 



where the system evolved hydrodynamically until hadronization and the 
hadron phase was described using a RQMD transport model. In this 
way the uncertainty of chemical non-equilibrium and its effects was cir- 
cumvented. The freedom in initial parametrization, however, was not 
explored, only the initialization sWN being used. Thus in the frame- 
work of a hydrodynamical description it is safe to say so far that the 
data seems to favour an equation of state with a phase transition, but 
no final conclusion can be drawn yet. 

From the hydrodynamical point of view elliptic flow at forward and 
backward rapidities is governed by the very same physics that creates 
elliptic flow at mid-rapidity. Thus to reproduce elliptic flow as a function 
of pseudorapidity all one has to do is to define how the initial density and 
shape of the system changes along the beam direction. This, however, 
would lead to a very counterintuitive initial shape of the system. The 
observed multiplicity stays approximatively constant within three units 
of pseudorapidity which - assuming that the longitudinal flow is ap- 
proximatively boost invariant - would mean that the initial entropy per 
unit of flow rapidity is approximatively constant within three units o f ra- 
pidity. On the other hand the elliptic flow data peaks at midrapidity ^4 
Combined with the requirement of constant entropy, this would mean 
that the initial system shape is most asymmetric at midrapidity and 
becomes more cylindrical when fragmentation regions are approached. 
It is difficult to imagine a physical process leadingto this kind of dis- 
tribution. So far hydrodynamical calculations ' — "I have not tried to 
find the initial shape to fit the data by trial and error, but have instead 
tried to formulate an intuitive parametrization for the initial shape and 
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calculated the elliptic flow based on that. These parametrizations either 
assume the shape not to change or to become more asymmetric towards 
the fragmentation regions. As a result the calculated elliptic flow fits the 
data only around midrapidity where \r)\ < 1 and stays almost constant 
in a large region of pseudorapidity. 

There is considerable freedom in choosing the initial state in hydro- 
dynamics and therefore it is premature to draw final conclusions as to 
whether the observed pseudorapidity dependence of elliptic flow really 
means a deviation from hydrodynamical behaviour. There may be a rea- 
sonable way to tune the initial conditions to reproduce the data which 
we have not thought about yet. 



5. Summary and outlook 

Hydrodynamical models have been surprisingly successful in explaining 
the flow data obtained at RHIC. The differential elliptic flow of charged 
particles is reproduced up to pt = 1-5 GeV, its centrality dependence 
works below impact parameter b « 6 fm and the differential elliptic flow 
of identified particles shows the mass ordering predicted by hydrody- 
namical calculations. That a model which is based on thermal averages 
is able to reproduce details at the 5-10% level is remarkable. This suc- 
cess, however, leaves us two complementary puzzles to solve: If the sys- 
tem does not thermalize, why is the observed elliptic flow so close to 
hydrodynamical elliptic flow? Or, if the system thermalizes, what is the 
mechanism responsible for it and how does one describe the thermaliza- 
tion process? At the time of this writing, there are no answers to either 
of these questions. 

The usefulness of flow studies and hydrodynamical models is not 
limited to deducing whether the collision system thermalized or not. As 
described in the previous chapter, flow is affected by the freeze-out tem- 
perature, the equation of state and the initial state. Thus it is possible 
to obtain information on all of them by studying flow, but, on the other 
hand, the complicated dependencies makes it challenging to constrain 
them individually. This requires careful simulations where one tries to 
reproduce as large an amount of the data as is possible and where the 
effects of all variables are taken into account. So far it is possible to say 
that the flow data favours an equation of state with a phase transition 
and initial state with densities well above the phase transition temper- 
ature. However, there are also open questions like the pseudorapidity 
dependence of flow and the HBT radii which must be solved before con- 
clusions can be drawn. 
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